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I Abstract. We consider non-rotating geodesic perfect fluid spacetimes which are 

purely radiative in the sense that the gravitational field satisfies the covariant 
transverse conditions divH = divE = 0. We show that when the shear tensor cr 
is degenerate, H, E and cr necessarily commute and hence the resulting spacetimes are 
hypersurface homogeneous of Bianchi class A (modulo some purely electric exceptions) . 
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1. Introduction 

(N \ 

^ . Although gravitational waves in perfect fluid spacetimes are usually studied via 
. transverse traceless tensor perturbations on a Friedman-Lemaitre-Robertson- Walker 
(FLRW) background, there is renewed interest in the so-called covariant approach [H El 
|3l H] in which gravitational radiation is described via the electric and magnetic parts of 
^ . the Weyl tensor, which represents the locally free part of the gravitational field. Perfect 
O ' fluid spacetimes in which the covariant transverse conditions divE = divH = hold at 
the non-perturbative level have been called purely radiative [5] and were shown to be 
^ . very restricted in the sense that they would have to obey two non-terminating chains 
of integrability conditions (see for example [HI [7]). All the evidence so far seems to 
indicate [8] that these conditions only hold in models with special symmetries, such as 
the spatially homogeneous spacetimes of Bianchi class A P,|TO]. In models with realistic 
inhomogeneity, the gravito-magnetic field would therefore necessarily be non-transverse 
at second and higher order. 

Recently we discussed purely radiative irrotational dust spacetimes for which the 
magnetic part of Weyl tensor is diagonal in the shear-electric eigenframe [H], as well 
as their generalization to purely radiative and geodesic perfect fluids with non-constant 
pressure [12]. In the latter paper, hereafter referred to as paper I, we showed that the 
Bianchi class A perfect fluids with non-constant pressure could be uniquely characterized 
(modulo some purely electric exceptions) as those purely radiative models for which the 
magnetic and electric part of the Weyl tensor and the shear tensor are simultaneously 
diagonalizable. For the case of constant pressure the same conclusion holds provided the 
spatial gradient of the matter density p vanishes: the velocity then becomes orthogonal 
to the p = constant hypersurfaces, and hence the flow is irrotational [T3l fl^ . Our aim 
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was to show that this analysis could be generalized to the case where H, E and cr 
are not necessarily simultaneously diagonalizable. The initial steps of this programme 
looked hopeful, in the sense that rational expressions could be constructed for all 
relevant quantities appearing in the description of the problem. However the sheer 
size of these expressions made any subsequent investigation of the ensuing integrability 
conditions virtually impossible. For technical reasons we therefore assume in the present 
contribution that the shear tensor is degenerate. A similar, but altogether different, 
assumption was made in |15|, where it was assumed that all kinematical variables, 
together with the electric and magnetic parts of the Weyl curvature were rotationally 
symmetric about a common spatial axis. This assumption gives rise to the so called 
partially locally rotationally symmetric (PLRS) cosmologies. The intersection of the 
PLRS cosmologies with our present class of models (i.e. those cosmologies having 
divH = divE = and a degenerate shear tensor), consists precisely of the LRS Bianchi 
cosmologies and the Kantowski-Sachs perfect fluids discussed in paper I. 

In §2 we give an overview of the results of paper I which remain applicable to the 
case under study, while in §3 and §4 we present a proof of the main result. 



2. Main equations 

We will follow the notations and conventions of paper I. All calculations were performed 
using the OFRAME package [IE] written in the symbolic computation language MAPLE 
and which is available from the authors. 

The energy momentum tensor of a perfect fluid has the form: 

Tab = (P + P)UaUb + pgab, UaU" = -1 (l) 

where p is the energy density in the rest-frame, p the pressure of the fluid and u the 
normalized 4-velocity. We also define E and H as the electric and magnetic part of the 
Weyl tensor with respect to the congruence 

Eab = Cacb,U^u' Hab = \eac,C^\eU^, (2) 

where Eabc = Vabcdu'^ is the spatial projection of the spacetime permutation tensor rjabcd- 
From the 'divH' Bianchi identity (see the appendix of [12]) and the assumption 
divH = we obtain [cr, E] = 0, implying the existence of a common cr and E eigenframe, 
which will be used henceforth. From the degeneracy of the shear tensor one can easily 
infer then that E must be degenerate too p^. Therefore we may assume that: 

cr = diag(-2a, cr, a) , . 

E = diag(-2E,E,E), ^ ^ 

where a and E are non-zero, as otherwise the space-time is FLRW. As explained in |12j . 
because all off-diagonal elements of E vanish, we have ^2 = fis = 0. By choosing the 
tetrad to be Fermi propagated we may also assume Qi = 0. Herewith the rotation rate 
of the tetrad about the ei-axis is fixed up to rotations 

^2 ©2 cos + 63 sin a 

63 —62 sin a -|- 63 cos a. 
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satisfying doa = 0. With ii = uj = the equations (4)-(12) of [12] simphfy to the 
following (Z is the spatial gradient of the expansion scalar, while n^.r^^qa are defined 
in terms of the rotation coefficients 7^^ = epsen^^ + d'^ap — dpas by Ua+i a-i = ('"a+5'«)/2, 
aa = {ra - qa)/'2 and n^a = n^+i + Ua-i): 

doO = -le^ - - ^{p + 3p) 

doP = -{P + P)0 (4) 



doa = a{cr-p)-E 



dia = -iZi + |(T(ri - gi) 



^2(7 = - Sag: 
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dsa = fZs + Scrrs 

d^E = l{ri-qi)E 

d2E = 3q2E 

d^E = -SrgE 



d,p = 

d2P = -9(^2 - 3aq2)a 

d^p = -9(^3 + Sars)^ 

-f^ii = 3(723 + ^2)0- 

H22 = -Sn^a 

= -3^2(7 

H12 = -(Zs + Sars) 

H23 = -f a(ri + qi) 

His = Z2- 3aq2 

while the commutators simplify to: 



(5) 



(6) 



(7) 



(8) 



(9) 



[do,di] = {2a-l9)di 

[do,d2] = -Ii3a + e)d2 

[80,83] = -li3a + e)ds (10) 

[^2, ^3] = {n2 + n3)di + ^3^2 + r2d3 

[83, 81] = rsdi + [rii + n3)d2 + qids 

[81,82] = q2di + ri82 + {ui + n2)83. 

Acting with [80, 81] on 6 and p implies qi = ri and Zi = 0, while evaluation of the 
following combination of commutators, 9a[82, 83](T + [82, c^sjp + 3a[82, ^3]^, leads to 

^2^3+^3^2 = 0. (11) 
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Furthermore, acting with [do, da] (a = 2, 3) on 9 gives us the time evolution of Z2 and 
Z3: 



doZ2 = -\aZ2 + f g2^' - ^^^^ 



J2 

d,Z; = -\aZ^ - f raa^ - QZ; 



which guarantees (cf. above) that we can rotate the (e2, 63) basis vectors such that 
for example Z2 = Z^ (as 9o(2'2/Z3) = 0). Note that Z2 and Z3 cannot both be 
zero, as then equation (fT2|) would imply q2 = = 0, such that by ([7]) the spatial 
gradient of p would vanish. The tetrad could then be further specified by rotating 
such that = ^{ri + qi) = 0, which by ([9]) would take us back to the case 
fcr, El = fcr, Hi = [E, Hi = |T2|. So henceforth we work in a fixed frame with 



Z2 = Z-^ ^ which by (fT^ implies r^ = —q2. Together with the Jacobi equations 
imply then 

dor, = {2a - h)r,, (13) 

doq2 = - ^2 + (5a - ^9)q2. (14) 

Applying now the commutators [9i,(92] and [93,9i] to 9 and combining the resulting 
equations, we obtain: 

2ri + 2ni +^2 + ^3 = 0. (15) 
Propagating Z2 = Z^ along 62 and 63 results in 

Z2ir2 + q-s) + ^a{nl - nl) = (16) 
which together with ffTSl) and the 'dotH' equations 

H<_ab> = cnilEab — 9Hab + 3ac<aHb^'^ 

leads to the key equation: 

{nl-nl)x = (17) 

where x = Z2{2E — Qa^) + 45g2a^. If we propagate x = twice along the fluid flow 
lines, we re-obtain Z2 = Z3 = 0. From (fT7|) we conclude therefore that n"^ — n"^ = and 
hence, by (fT6l) . q^ = —r2. The resulting cases n2 + = and ^2 + n3 7^ (and hence 
^^2 = '^s 7^ 0) will be investigated in the following paragraphs. 

3. The case n2 + ns = 

The 'dotE' equations 

E<:ab> = CUrlHab - 9Eab + 3(yc<aEb>'' - \iP + P)^ab 

imply that also curlH is diagonal in the shear eigenframe. Using the off-diagonal 
components of curlH we can propagate n2 -|- ns = along 62 and 63 respectively, 
which leads to 

{Z2 - Aq2a){n2 - n) = = {Z2 - Aq2a){n2 + n). (18) 
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First notice that Z2 — 4^20" = would lead to a vanishing spatial gradient of the 
matter density. In fact, propagating Z2 — ^q2<y = along the fluid flow lines gives 
q2{8E — Tcr"^) = 0, after which the propagation of 8E — la'^ = along 62 would 
yield 52^"^ = and hence Z2 = = 0. From (fTSl) we conclude therefore that 
n2 = ris = vi = 0. 

Propagating n2 = along ei one obtains then 

E-|-2a2- — + - + 2r2g2 = 0. (19) 



Acting with the commutator [(?o,(?2] on Z2 and using (fT9|) allows us to eliminate E, 
which implies 

- 12r2g2' + (2p + + 2ea - h')ql - + ^Z^ = (20) 

and which also shows that ^2 7^ 0. Propagating the latter relation along 63 and using 
(1201) to eliminate p one finds: 

{Z2-3q2a){4Z2-21q2a) =0. 

When we propagate the first factor along the fluid flow lines, we obtain an expression 
for p which, after substitution in ( |T9|) . would lead to the "anti- Newtonian" result E = 0, 
which is impossible . On the other hand, if we apply the same procedure to the second 
factor of ([3]), we get (j(4Z2 — 9g20") = which gives a contradiction with 4^2 — 21g20" = 0. 

4. The case n2 + ns 7^ 

We now have = n2 ^ 0, which, when propagated along ei, implies 

E = 2q^q2 - 2r^n2 + \p + 2a^ + lad - -n\- r\. (21) 

Acting with the commutator [So, d'2\ on p, Z2 and r\ we arrive at the following expressions 
involving the matter density and the pressure: 

- yg20-^2 + ISnariO-^ + 4n2rip - ^n2ri9^ + 5Z^ + 

2 

+ 12g^g3 - -e^qj - I2n2riql + 2pg^ - Qnlql - Grjql + An2riae 
o 

+24n2rig3g2 + 2qlae - V2rxn\ - 2Ar?^r\ - \2n2r\ + 45(tV^ = 



(22) 



30Z2n2p - 18riZ2p + 54r^Z2 - 90Z2n^ - 18riZ2cr0 - 108riZ2g3g2 
+30^2^2^^ + 180^2712^3^2 + 2025a^g2^2 + 1377rig2fT^ - 405riZ2a^ 
+6riZ2^^ - 126riZ2n^ + \8r\Z2n2 - 333^2^2^^ - 10^2^2^^ = 



(23) 
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and 



27 

p = - - 6a0 - 4p - 18g3g2 + + 9nl + ISri^a + 9rj 
189q2a^ + 54q2a^e + 36g2(xp + 216g|(Tg3 - 12g2(T^2 

^ 2Z^ 

108q2crnl + 216g20"ri'n-2 + 108^20^^1 



(24) 



2^2 

Eliminating ^3 from fl22]) and fl25]) gives the relation 

- 50^1^2 + 918rig^cj^ + 30Z|ri + 1350cT^g^n2 + 1836r^g2a^n2 
+2700(7^^2^2^1 + 153riZ2(T2g^ _ goOrlZ2Cx^n2 - 921Z2n2(y'^ql 
-2mq2(TZlri + 34552^x^2^^2 + 270^^5*^2 - 594Z2n^aVi = 0. 

(25) 

Notice that ri has to be different from zero (see the appendix), while (1141) guarantees 
that q2 7^ 0. This enables us to introduce new dimensionless variables, 

Z2 q2 , n2 n qs P 

X = , y = — , n = —,n = — , r = — ,q= — and m = — , 

q2(7 ri a ri a a 

which satisfy the following evolution equations 

a~^doX = - [{n + lfr^ -2qyr + {2h^ -6h-6m + 153)/18]x + 4:5/2 

a^^don = —6n 

a-^dor = -{n + lfr^ + 2qyr^ + -{h^ -6h-3m-27)/9r 

(^'^doy = -{x~3) y 

and 

a^^d2X = 2r[—^yx'^ + {n + l + 2y'^)xy^^—6ny^^ — 3y] 

a^^d2n = ry [(I — I n — I n^) X + 4n^ — 2] + 2 gn 

a^^d2r = (5 — 4n — I X + I xn) r^y — 2 gr 

cr^^d2y = (^^ X — ^ xn + An — 3) y'^r + 3 qy — 2rn — 2r. 

Hereby flM|) was used in order to substitute for all occurrences of the pressure in the 
right hand sides. In what follows we will show that the propagation of equations ( !23l) 
and (125!) implies that x is constant. Then however (p7|) (a) and dQd2X = c?g92X = 
simplify to {x — 6){x — 3)^y^ + 9x = and a{x — 6){x — 3)^?/^ = 0. This implies x = 0, 
which brings us back to the degenerate shear case of paper / and thereby ends the proof 
of our theorem. 

First we rewrite equations ( 123|1 and ( l25l) as follows 

— 18 mx + 30 mxn — 18 /ix — 108 qyrx + 30 hxn + 180 xyrnq 
+54 r\ - 90 xr^n^ + 2025 n + 1377 - 405 x + 6 h^x - 126 r^n^x 
+18 r^nx — 333 nx — 10 xnh'^ = 

(28) 



(27) 
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1350 y^n + 918 y^ - 50 x^y^n + 30 x^y^ + 1836 n + 2700 

+ 153y^x - 927xy^n - 207 ?/V + 345?/Vn + 270x 

-gOOnx - 594xn^ = 0. (29) 

Propagation along 62 of ( !28l) and ( !29|) results in a linear homogeneous system in q 
and r, the coefficients of which are polynomials in n, x and y. Eliminating n from the 
determinant of this system and from (l29l) yields a polynomial relation in x and y: 

Pi(x,?/) = -545500000 x^^y^^ + 22883025000 x^^y^^ - 12001000000 x^^y^^ 
+ ... (76 terms) . . . + 7746297620889600 x^ + 140466196858798080 y^ 
-20656793655705600 X = 

(30) 

We construct a second polynomial relation between x and y by propagating equation 
( l29l) along the fluid flow, which yields 

- 300 nx^y^ + 250 n^x^y^ + 90 x^y^ - 3600 n'^x^y^ + 8460 nx^y^ 
-3780x^1/^ + 11925nV?/2 - 54900 nxV + 28647 xV 
+ 144018 nx?/2 - 69714 x?/2 - 8748 n^xy^ + 2718 nV 
+2970 n V + 810 x^ - 4050 nx^ + 38556 y^ - 161028 ni/^ 
- 12960 X - 49896 ri^x + 1296 n^x + 65880 nx + 162000 
-127656 - 162648 n + 38880 = 

(31) 

Eliminating n from ( |29l) and (1311) results then in 

P2(a;, y) = 10721500 x^ V - 380191500 x V + 6061903695 x V 
+ ... (27 terms) ... - 1402072053120 x - 1800380390400 = 0. 

(32) 

One now can calculate the resultant of (130|) and (132|) with respect to y: this yields a 
polynomial in x which is not identically and thereby shows that x is a constant. 

5. Conclusion 

A geodesic and non-rotating perfect fluid (where 'non-rotating' can be dropped when we 
have non-constant pressure), with degenerate shear tensor cr and with divE = divH = 
has commuting H, E and cr. By paper I this implies that the resulting spacetimes are 
spatially homogeneous of Bianchi class A. More particularly, because of the degeneracy 
of cr, they are of type VIq (in the non-LRS case), or, in the LRS case, of Bianchi 
types I (VIIq), II, VIII or IX. The only exceptions arise in the purely electric case, 
where also the pseudo-spherically symmetric Bianchi class B type III spacetimes and 
the Kantowski-Sachs spacetimes are allowed. 
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6. Appendix 

the case ri = 

Propagation of ri = along 62 gives n2(24g20" — 5Z2) = 0. As the evolution of 
24^20" — 5Z2 = along the fluid flow would lead to a = we have n2 = 0. Propagating 
( I22|) along 62 and substituting n2 = gives {Z2 — 3g2cr)(4^2 — 21^20") = 0. The first 
factor must be nonzero, as otherwise (with ri = 0) we would have a diagonal H. When 
the second factor is we obtain from ( |22l) 

AO [X 2 

— or^ + 12q;q2 - -e^ + 2p + 2ae = 0. (33) 
lb 3 

Eliminating p from the time evolution of 4^2 — 21g2C = and (1331) gives a = 0. 
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